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ABSTRACT 


A finite element formulation was used to solve the steady- 
state two-dimensional conduction heat transfer equation in 
the condenser wall section of an internally finned rotating 
heat pipe. A FORTRAN program using this method was coupled 
with the COPES/CONMIN program for sensitivity analysis of 
design variables and for automated design of the internal 
heat pipe geometry. 

With water as the working fluid, numerical results ob- 
tained for copper and stainless steel heat pipe condenser 
sections indicated that for the maximum heat transfer rate, 
the designer should machine as many fins as the condenser 
material and the manufacturing process will allow. A saw 


tooth profile is preferable to spacing between fins. 





ir. 


ei. 


ei. 


ney 


Waris Of CONTENTS 


INTRODUCTION- - = - - = -=-2+efee-77-772722e2+0¢- 14 
A. THE ROTATING HEAT PIPE- - - - -------- - 14 
EeOeErATING EBIMETS OF A ROTATING HEAT PIPE : 
CONDENSER SECTION - - - - ---+-+e------- is 
1. The Sonic Limit - ----+-+--+------ iS 
2. Boiling Limit - ----+--------2-- Hue 
3. Entrainment Limit - - ----+------e-- ig 
4. Condensing Limit- - ------+-+e----- alee 
C. ANALYSIS OF AN INTERNALLY FINNED 
ROTATING HEAT PIPE- - - - - - - ---+-+--e-e-- 20 
D. THESIS OBJECTIVES - - ---*+-+e4-e*-+-+e-+-+-- 21 
NUMERICAL OPTIMIZATION- - - - ----+-+-+-+--- - 29 
A. BACKGROUND- - - -----+-+-+-----+---- 22 
B. CONMIN- - - - - ----+-ee-e--e7------ - me 
M@eeGOPES - = == ------------+------ 24 
MP Gee Iiimic GUEDELINES= - - --+----+.44- 24 
FINITE ELEMENT SOLUTION - - - - ----+-+e-+-+-- - 28 
A. REVIEW OF PREVIOUS ANALYSIS - - - ----+-+- - Ze 
B. COMPUTER PROGRAM DEVELOPMENT- - - - - ---- - 30 
C. SENSITIVITY ANALYSIS- - - ----+-+-+----- 32 
D. DESIGN OPTIMIZATION - - ----+-+--+-+--- - S 
RESULTS - - ----+--+-+-+e------4--42+-2-- au 
 ROMUCMIO N= == = - = = = - - - - = - + - ou) 
B. SENSITIVITY ANALYSIS- - - ----+-+-+-+-+--- ay 
C. CONSTRAINED OPTIMIZATION- - - - - ------ - So 
D. UNCONSTRAINED OPTIMIZATION- - - - ------ - LO 





E. A CAUTIONARY NOTE - - - - ---+-+e-+-+---- - 43 


V. CONCLUSIONS - -----+-+-+e---+-7-+-7+4-+-+-+-2-- 4S 

VI. RECOMMENDATIONS - - ----+-e-+e-+-+e-+---+---2e-- 47 
VII. FIGURES - ----+-+e-+e-ee-e-4-4---4-4-4-4-4--%2-- 4 8 
APPENDIX A - MANUAL FOR USING COMPUTER PROGRAM- - - - - - 60 
APPENDIX B - PROGRAM LISTING- - - - ----+-+-+-+--- - G1 
BIBLIOGRAPHY- - ----+-+-+e-+-e4---4-4-4-4-4-4-42-2-- eae 
INITIAL DISTRIBUTION LIST - ----+-+--+-+--+---- - 106 





ket OF TABLES 


Specification of a Typical Rotating Heat Pipe- - 


Basic Design for Sensitivity Analysis- - - - - - 


Constrained Optimization Results for 
Stainless Steel, 


Copper, 


and Ceramic Material- - 


Unconstrained Optimization Results for Copper, 


Stainless Steel, and Ceramic Material- - - - - - 


Comparison of Unconstrained and Constrained 
Optimum Design for a Copper Condenser Section 
with External Heat Transfer Coefficient of 


20000 


BTU/HR : 


rae 


ca 


18 


a0 


od 


Eo 


43 





im). 


ial. 


ie 2 


A-l 


fit. OF PaGuRes 
Schematic Drawing of a Rotating Heat Pipe- - - -- - 48 


Operating Limits of a Typical, Water- 
Filled Rotating Heat Pipe- -----+--+--+-- - Lg 


Internally Finned Condenser Geometry, 
Showing Fins, Troughs, and Lines of Symmetry - - - - 50 


Condenser Geometry Considered with 25 
Linear Triangular Finite Elements- - -------- a 


Differential Equation and Boundary Conditions 
Considered in the Analysis of Purnomo [l]J- - - - - - 52 


Comparison of Trough Condensate Film 


Thickness vs. Length Along the Condenser - - - - - - a3 
Heat Transfer Rate (Q) vs. Number of Fins- - - - - - 54 
Heat Transfer Rate (Q) vs. RPM---+-+e--+e-+--- - 55 
Heat Transfer Rate (Q) vs. External 

Heat Transfer Coefficient- - -----+-+----- - 56 
Heat Transfer Rate (Q) vs. Condenser Length- - - - - oi 
Heat Transfer Rate (Q) vs. Cone Half Angle - - - - - 58 
Heat Transfer Rate (Q) vs. Fin Half Angle- - - - - - 59 


Specification for Input Data to Determine 
the Coordinates of the System Nodal Points - - - - - 90 





A 


A 
S 


b (BVIN) 


ee (HFG) 


Ke Cr ) 


h 


ia (CW) 
W 


L 


M (AMTOT) 
N 

ie 
aecQTOT) 
Q 


S 


R (RBASE) 


Rabie Or SYMBOLS 
cross sectional area for flow in aoe finite 
element area 
2 


iicsiGemcimadee agreaq of a2 Smooth tube in ft 
Morten Enomrineine his Linite element factor 
sonic speed in ft/sec; finite element factor 

Z 


acceleration of gravity in ft/hr 


convective heat transfer coefficient in 
Btu/hr-ft -°F 


latent heat of vaporization in Btu/lbm 


thermal conductivity of condensate film in 


Btu/hr-ft-oF 


thermal conductivity of condenser wall in 


Btu/nr-ft-TF 

finite element sides 

mass flow rate of condensate in lbm/hr 
two-dimensional linear shape function 
pressure of the vapor in lbf/f£t? 

heat transfer rate in Btu/hr 


heat transfer rate through a smooth tube in 
Bru/ hr 


internal radius of condenser in ft; thermal 
resistance in hr-f/Btu 


temperature 


Vieloe@iy on) liquid in ft/sec 
axis of Cartesian system coordinate 


coordinate measuring distance along the 
condenser length 





SA 

2 

Shik K 

a (FANGL) 
6* (DEL) 
mC EPS ) 

fee PHT) 
Ps (RHOF ) 
Py 

Oo 

ee Ur) 
Py 

w (OMEGA) 


coordinate measuring distance perpendicular to fin 
Ssuriace 


coordinate measuring distance along fin surface 


fin half angle in degrees 

econdemsece t1im thickness in ft 
MOcatmeroueh WHat Li ft 

condenser cone half angle in degrees 


denctnweonr=tie liquid. in ibm/ ft? 


density of the vapor in lbm/ft> 
Uma @eommmenctoneor the I1iguid in 1bf/ft 
ViseOsimagot themiiquad in 2lbm/ft-sec 


viscosity of the vapor in lbm/ft-sec 


angular velocity rad/sec 


EO 





ALFA(a) 
BFIN(b) 


BVIN 

CALFA 
CANGL 
SBASE 
SEXIT 


er (X,) 


CL 
eit 
ero 
eit 
DIV 
mMrOT 
DOBF 
mor H 
FANGL 


HFG 
aE |, 
JLC 


aC 


KFF 


ir LN 


KT 


TABLE OF FORTRAN VARIABLE NAMES 


FIN HALF ANGLE (RADIANS) 

Hemet OF BIWf CINCHES ) 

EGE OF@rIN (FEET) 

COSINE OF ALFA 

GCONERHALE@ANGEL (DEGREES) 

iS vero tneOrernRarer OF CONDENSER (FEET) 

Pfeil iE CIR@UMPERENCE @Al CONDENSER EXIT (FEET) 


Pen eeecONDUCTIViTy OF CONDENSATE FILM 
CorU/AR Else. F) 


CONDENSER LENGTH (FEET) 

CONDENSER LENGTH (INCHES) 

Ceri Of FHI 

Col VERGENGCE CRITERION 

FLOATING POINT VALUE OF NDIV 
CONDENSATE MASS FLOW RATE 

NUMBER SOF COLUMN WITHIN THE FIN 
NUMBER OL SCOLUMN WITHIN THE TROUGH 
FIN HALF ANGLE (DEGREES) 


CCW ECE LY & HEAT De meee Orr: LOLENT 
CoaUyannewn l= <° F) 


LATENT HEAT OF VAPORIZATION (CBTU/LBM) 
PoE ween T NUMBER 


NMUMBeweOb ssYotEM NODAL POINT LOCATED AT 
eeewcrn ben. OF) SY SlEM COORDINATE 


Nek ete srothM NODAL POINT LOCATED AT 
fiero OF THE SYMMETRY BOUNDARY 
ites ObetNTERSECTION BETWEEN THE 
FIN AND THE CONDENSER WALL 


NOE Or  SovorEM NODAL POINTS LOCATED ALONG 
ieee te CONVECTIVE BOUNDARY 


NOSE Reba colEli NODAL POINTS LOCATED ON 
THE SYMMEFRIC BOUNDARY OF TRIANGULAR FIN 
SECTION NOT COUNTING POINTS AT BASE AND APEX 


NUMBER OF ROWS WITHIN THE WALL SECTION 


el 





NBAN 
NBOTF 
NBOTI 
NDIV 
NEFB 
NEL 
feot 
NSNP 
PHI($) 

- Ge) 
RBASE 
RBASEI(R) 
REXIT 
SALFA 
SseHt 

ia iCK 
ea iCK LT 
err 
erin 
ZOAC c/a) 


SYolEM BAND WIDTH 

Pot been Abe BOTTOM SIDE 

Pavel ELEMENT Ar BOTTOM SIDE 

NUMBER SOF INCREMENT 

ELEMENT NUMBER AT BASE FIN 

(MBean OF ELEMENTS 

ELEMENT NUMBER AT END OF TROUGH 

UME OF eee l EM NODAL POINTS 

CONE HALF ANGLE (RADIANS) 

El 

Mis eeaw ole ot CONDENSER BASE (FEET) 
iS oes RAD mUSe OF GONDENSER BASE (INCHES) 
ieee eis Or CONDENSER EXIT (FEET) 
SNE OF ALTA 

SNE OF Tene 

CONES Eke WAL THICKNESS (FEET) 

CO MENSReeWAli THICKNESS CINCHES) 
TANGENT OF PH 

NUMBER OF FINS 

RATIO OF TROUGH WIDTH TO FIN BASE WIDTH 


tz 





ACKNOWLEDGMENTS 


The author wishes to express his appreciation to Dr. P. 
OW. Marto, Professor of Mechanical Engineering, Dr. D. Salinas, 
Associate Professor of Mechanical Engineering, and Dr. G. N. 
Vanderplaats, Associate Professor of Mechanical Engineering 
for their advice and guidance throughout the development of 
this thesis. 

Finally, the author wishes to thank his wife, Mary, for 
her understanding and encouragement during the course of this 


eeeudy. 


i 





Pee ODUCT LON 


fee lH ROTATING HEAT PIPE 

The rotating heat pipe is a closed container designed 
to transfer a large amount of heat in rotating machinery. 
Its three main component parts are: a cylindrical evapora- 
Or, a truncated cone condenser, and a working fluid as 
shown in Figure l. 

At rotation above the critical speed of a rotating heat 
Dipe, the working fluid forms an annulus in the evaporator, 
and will be vaporized by heat addition to it. The vapor 
flows toward the condenser as a result of a pressure differ- 
ence, transporting the latent heat of vaporization with it. 
External cooling of the condenser causes the vapor to con- 
dense on the inner wall and release its latent heat of eva- 
poration. The centrifugal force due to the rotation has a 
component acting along the condenser wall that will act to 
drive the condensate back to the evaporator where the cycle 
ls repeated. 

In a conventional heat pipe, the force driving the con- 
censate back to the evaporator is due to capillary action, 
which poses a limit to its operation. The rotating heat 
Pipe is not limited by capillary action and, unlike the 
thermosyphon which depends on gravity to cause condensate 


return, can be used in any orientation [1]. 


14 





See OPERATING LIMITS OF A ROTATING HEAT PIPE 

The first theoretical investigation of the rotating 
heat pipe conducted at the Naval Postgraduate School was per- 
formed by Ballback [2] in 1969. He studied the limitations 
of performance imposed on the rotating heat pipe due to 
various fluid dynamic mechanisms. Using existing theory and 
experimental correlations, he was able to estimate the sonic 
[emit, boiling limit, entrainment limit, and the condensing 
limit of performance. 

Mee the Sonic Limit 

When increasing the heat flux in a rotating heat 

pipe, it is possible to reach a limiting flow rate of the 
vapor brought on by a choked flow condition in the pipe. 
mes Condition imposes a limiting value on the amount of 
energy the vapor can transport, thus reducing the effective- 
ness of the heat pipe. The limiting heat transfer rate 


becomes 
pe Opty She, Gl 


and the vapor velocity is considered to be sonic, 


ie - co = y BkRT G2) 
where 

a wevelLOcTuy Of the vapor in ft/sec, and 

A = cross sectional area for the vapor flow in ft’ 
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Ee pecOmmem7eClOCzEyY in ft/sec 
Ey = Gha weed lOnal constant, 32.1739 ft-lbm/1bf-sec* 
k = patio of specific heats 
R = gas constant in ft-lbf/lbm °R and 
z = absolute temperature in OR. 
fe SOLling Limit 
Kutateladze [3] postulated that the transition from 


hucleate to film boiling is totally a hydrodynamic process. 


He determined a theoretical formula for predicting the burn- 


@ae flux 
Mee- K 7p A. h._{o e(p- - oe ee C2) 
it ioe oo f 4 
where 
K = constant value 
= density of the vapor in nee 
A, = heat transfer area in the boiler in ft< 
Ne, = latent heat of vaporization in Btu/lbm 
Oo = surface tension in l1bf/ft 
g = aeceleration Of gravity in ft/hr? 
Pe = @demsity or 2 lturd in Lbm/ ft? 
- = density of vapor in Lbm/ ft? 


The experimental data obtained by Kutateladze suggested a 


value for K in the range of 0.13 to 0.19. 


nee 





fe EMtreainment Limit 
The flooding constraint in a wickless heat pipe was 
examined by Sakhuja [4] who developed the correlation 
ACh. Vv eDCp =p.) 0. 
Q, = x ae To eas a 
ilar Co soe) }? 
where 
Q, heat transfer rate in Btu/hr 
ao flow area in ft? 
C dimensionless constant, 0.725 for tube with sharp 
edged flange 
Deg latent heat of vaporization in Btu/lbm 
g acceleration due to gravity in ft/hr- 
D inside diameter of heat pipe in ft 
ai density of the fluid in lbm/ft? 
oe, density of the vapor in Lbm/ ft? 


4. Condensing Limit 


Pauuback (2) determined the condensation solution 


for a rotating heat pipe by modeling the condenser section 


of a rotating heat pipe as a rotating truncated cone. He 


developed the following expression for the condensation 


Seirmat : 


i 
2 3 = 8 8 3 
jae Betis yt .* aoe 
1t{= - } ee + sin o J -R ~} 
S oe O O 
Us Sin @¢ 


co 


coy 





EOmoteat tEranster rate in Btu/hr 


GiermalwecOnGauectcivity Or the condensate film 
imepeuyne—ft-OF 


density of fluid in Lbm/ ft? 

apewlar velocity imei/hr 

Patent Neat of Vapomazation in Btu/lbm 
saturation temperature in °F 

Mcrdeswall temperature, an, 2 F 

Viseoowty Of flutd insibm/ft—-hr 

half cone angle in degrees 

minimum wall radius in ft 

length along the wall of the condenser in ft 


viscosity of the fluid in lbm/ft-sec 


The condensing limit equation is a function of the geometry 


and speed of the rotating heat pipe, and the physical pro- 


Mertites of the working fluid. 


Seiieeaicil §5)] calculated these limitations for a heat 


Pipe with specific physical characteristics as shown in 


Table I, with the results shown in Figure 2. 


done, 


BeCciLiGation Of 4a Lypical Rotating Heat Pipe 


Length 14.000 aLigvelloveste! 
Minimum diameter 2.008 inches 
Wall thickness Ot, dees inches 
Mmserndimahalf angle OG degree 
Rotating speed ZING RPM 
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Q@bviously from the results in Figure 2, the condensing limit 
is the predominant limitation for the amount of heat that can 
be transferred from the heat pipe. However, the other 
limitations may become important as the heat pipe geometry and 
operating conditions are varied. 

In order to augment the heat transfer capacity of the 
heat pipe, recent efforts have been aimed at raising the con- 
densing limit line which may be accomplished by: 

@. a high value of cone angle, to increase the centrifugal 
eruving force’, 

b. some type of promoter of dropwise condensation to in- 
crease the value of the inside heat transfer coeffi- 
egemt, h, or 

ec. use of an internally finned condenser to increase 
the inner wall surface area and the value of h, since 
the presence of a fin will decrease the effective 
condensate film thickness. 

A high value of cone angle means a departure from the 
cylinder or shaft shape. Since the principal known appli- 
@ati0n for the rotating heat pipe is in the cooling of rota- 
ting machinery this approach was not pursued. Although 
effective promoters of dropwise condensation exist, none as 
yet can be considered to be permanent, and this approach was 
likewise ruled out. The remaining alternative, using inter- 
nal fins in the condenser section to raise the condensing 


limit, was seen as the best choice. 


IL 





PeerNNALYSiS OF THE INTERNALLY FINNED ROTATING HEAT PIPE 
Pursuing the addition of internal fins as a way to raise 
the condensing limit, Schafer [6] developed an analytical 
model for a heat pipe with a triangular fin profile as shown 
in Figure 3. He assumed one-dimensional heat conduction 
through the wall and fin. Corley [7] for this same case 
developed a two-dimensional heat conduction model using a 
Finite Element Method, and also assumed a parabolic tempera- 
ture distribution along the fin surface. His results indi- 
cated a significant improvement in heat transfer performance 
of about 75% above that predicted by the one-dimensional model 
of Schafer [6]. However, Corley [7] cautiously noted a 
probable error of 50% existed at the fin apex, and conse- 
quently mentioned that there may be a total heat transfer 
error of as high as 15%. Tantrakul [5] modified Corley's 
computer program by increasing the number of finite elements 
in order to minimize the heat transfer error at the apex of 
the fin. His results with this modification converged with 
the results of Corley. Purnomo [1] developed a two-dimen- 
Sional Finite Element Method solution using a linear trian- 
gular finite element model as shown in Figure 4. Purnomo's 
[1] Finite Element Method program also worked and converged. 
Purnomo's [1] code, when made to approach the geometry of a 
smooth tube, did not agree with the analytical and experi- 
mental data obtained by Schafer [6] for a smooth tube. This 


Cast doubts about the validity of Purnomo's cade. 
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The parametric studies conducted using Purnomo's code 
gave no clear indication of the best condenser geometry to 
Meaeimize heat transfer. Also, his code was tedious to use 
and required numerous runs to obtain data since it was 
written to perform only one analysis at a time. 

Clearly a computer program that could make numerous runs 
with minimal data input and could also automatically find 


improved designs would be valuable. 


See laeSslS OBJECTIVES 
The objectives of this thesis were therefore: 


1. To modify Purnomo's [1] computer program so that it 
is compatible with the COPES/COMNIN program [8] and 
can be used for analysis and automated design of 
rotating heat pipes (internally finned or smooth). 


2. To compare results using Purnomo's code with analyti- 
cal results for a smooth tube obtained by Schafer [6] 
to determine if and where an error exists. 


eeelo use the sensitivity analysis capability of COPES/ 
CONMENSin conjunction with the modified program to 
study heat transfer in an internally finned rotating 
heat pipe. 


4. To use the resulting program to obtain an optimum 
design for an internally finned rotating heat pipe 
to obtain experimental data to compare with the 
analytical results. 


5. To use the resulting program to obtain numerical 


results in place of data obtained from expensive 
experimental operations. 


Zale 





ieee Witenes eee LMT ZATION 


A. BACKGROUND 

Most design processes require the minimization oY maxi- 
mization of some parameter which may be called the design 
objective. For the design to be acceptable, it must satisfy 
a variety of physical, aesthetic, economic and, on occasion, 
political limitations which are referred to here as design 
constraints. While part of the design problem may not be 
easily quantified, most of the design criteria can be 
described in numerical terms. 

To the extent that the problem can be stated in numerical 
terms, a computer program can be written to perform the 
necessary calculations. For this reason, computer analysis 
is commonplace in most engineering organizations. For 
example, in structural design the configuration, materials, 
and loads may be defined and a finite element analysis com- 
puter code is used to calculate stresses, deflections, and 
other response quantities of interest. If any of these 
parameters are not within prescribed bounds, the engineer 
may change the structural member sizes and rerun the program. 
The computer code therefore provides only the analysis of a 
proposed design, with the engineer making the actual design 
decisions. This approach to design, which may be called 


computer-aided design, is commonly used today. 
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Another common use of analysis codes 1s in tradeoff 
studies. For example, an aircraft trajectory analysis code 
may be run repetitively for several payloads, calculating 
the aircraft range, to determine the range-payload sensitivity. 

A logical extension to computer-aided design is fully 
automated design, where the computer makes the actual design 
decisions, or performs trade-off studies with a minimum of 


mean=machine interaction [9]. 


B. CONSTRAINED FUNCTION MINIMIZATION (CCONMIN) 

Vanderplaats [10] developed an optimization program CONMIN 
capable of optimizing a very wide class of engineering pro- 
blems. CONMIN is a FORTRAN program, in subroutine form, that 
optimizes a multi-variable function subject to a set of in- 
equality constraints based on Zoutendijk's [11] method of 
meesdble directions [12]. 

Three basic definitions are required to discuss the 
Mee Of CONMIN: 

Design Variables - Those parameters which the optimization 
program 1s permitted to change in order to improve the 
design. Design variables appear only on the right hand 
Side on an equation and are continuous. 

Design Constraints - Any parameter which must not exceed 
specified bounds for the design to be acceptable. 

Design constraints may be linear or nonlinear, implicit 


Peeexplicit, but they must be continuous functions of 


De 





the design variables. Design constraints appear only 
on the left side of the equations. 
Objective Function - The parameter which is going to be 

minimized or maximized during the optimization process. 

The objective function may be linear or nonlinear, 

Pict Oreexplicic, ama must be a continuous function 

of the design variables. The objective function usually 

appears on the left side of an equation, but it may 

appear on the right side if it is also a design variable. 
Design constraints and objective functions are usually inter- 


changeable. 


Bec ONLTROL PROGRAM FOR ENGINEERING SYNTHESIS (COPES) 

Recall that the optimization program, CONMIN, was written 
in subroutine form. Vanderplaats [8] has developed a main 
program to simplify the use of CONMIN and to further aid in 
the optimization process. The user must supply an analysis 
subroutine named ANALIZ. What follows are programming guide- 


lines to ensure compatability with COPES. 


D. PROGRAMMING GUIDELINES 

In developing any computer code for engineering analysis, 
mms prudent to write the code in such a way that it is 
easily coupled to a general synthesis program such as COPES. 
Therefore, a general programming practice is outlined here 


which in no way inhibits the use of the computer program in 
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Meer traditional role as an analytic tool, but allows for 


memple adaption to COPES. This approach is considered good 


meogramming practice and provides considerable flexibility 


of design options. Only five basic rules must be followed: 


ft. 


Te 


Write the code in subroutine form with the primary 
routine called as SUBROUTINE ANALIZ(ICALC). The 

name ANALIZ is compatable with the COPES program and 
mesWe ers cl caoleulartion GomtrolseeNote that subroutine 
ANALIZ may call numerous other subroutines as required 
to perform the necessary calculations. 

Seomene the procram into INPUT, EXECUTION, and OUTPUT. 
five, calculationwcontrol., LECALC, will determine the 
portion of the analysis code to be executed. ICALC=1; 
the program reads all data required to perform the 

oes cee so, any initialization of constants 

which will be used repetively during execution is 

done here. This initial input information is printed 
here for later reference and for program debugging. 
ICALC=2; the program performs the execution phase of 
the analysis task. No data reading or printing is 
done here, except on user-defined scratch dise. Data 
may be printed here during program debugging, in 
which case it should be controlled by a print con- 
trol parameter which is read during input. [In this 


way, this print may be turned off after the program 


ZS 





eT. 


ny . 


is debugged, but may be used again during future 
program expansion debugging. The reason that print- 
ing is not allowed during execution is that when 
optimization is being done, the code will be called 
many times with ICALC=2, resulting in voluminous 
print. ICALC=3; the results of the anlaysis are 
printed. Also the essential input parameters which 
may have been changed during optimization should be 
printed here for easy reference. In summary, when: 
ICALC = 1 Read input data. 

HEALC ="2 Execute the analysis. 

Gene = so Prant the results. 

Store all parameters which may be design variables, 
objective functions or constraints in a single 
labeled common block called GLOBCM. The order in 
which they are stored 1s arbitrary. A listing of the 
COPES program should be checked to see how many para- 
meters may be stored in GLOBCM (the dimension of 

ee ee ee cise seribution Of COPES allows for 1500 
parameters. 

During execution or output, no parameters which are 
read during input should be updated. For example, 

if variable X is initialized during input, the 
execution segment must not update X such as X=X + 
Gees instead a new variable, Y=X + 3.2 must be 


@erined. 
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V. Write all programs in standard language, avoiding 
machine dependent capabilities such as seven tetter 
Beer eenemnes Coben wWhtle tmrs guideline is not 
essential to the use of the analysis code within the 
COPES program, it makes the analysis code much more 
transportable between different computer systems, a 
capability which easily justifies a slight reduction 
in efficiency on a given machine. 

Adherence to these guidelines not only leads to a more 
readable and machine independent computer code, but allows 
mos Code to be coupled to the COPES program without modifi- 
@aetlon. 

Having written the analysis code, it may be executed 
elther with a simple main program or within the COPES pro- 
gram to perform the analysis. To insure that guideline IV 
is followed, the following main test program is recommended. 
Note that this program calls ANALIZ twice with ICALC=2 and 
ICALC=3, to show that the same result is obtained repeti- 
ively. 

C MAIN PROGRAM TQ CHECK SUBROUTINE ANALIZ. 
c Pee, EXECUTE, AND PRINT 
berto ICALC=1,3 
CALL ANALIZ (CICALC) 
EeecULE AND PRENT AGAIN TO BE SURE THE RESULTS 
DO NOT CHANGE 
we 20 ILCALC=2,3 
20 CALL ANALIZ (CICALC) 


SLOP 
END 


CHO 
OQ 
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ii eke Cuewent SOLUTION 


See VIEW OF THE PREVIOUS ANALYSIS 


Schafer [6] studied the one-dimensional model heat trans- 


fer solution and Corley [7] studied the two-dimensional 


model for an internally finned rotating heat pipe. Both 


used the same assumptions and boundary conditions based upon 


the analysis of Ballback [2], which are similar to those used 


in the Nusselt analysis of film condensation on a flat wall. 


The more important of those assumptions are: 


i. 


ee 


steady state operation, 

film condensation, as opposed to dropwise condensation, 
laminar flow of the condensate film along both the 
fin and the trough, 

Seatre Dalance of forces within the condensate, 
one-dimensional conduction heat transfer through 
the film thickness (no convective heat transfer in 
the condensate film), 

no liquid - vapor interfacial shear forces, 

no condensate subcooling, 

zero heat flux boundary conditions on both sides of 
the wall section (symmetry conditions), as shown 

fay PLOUre o> 


saturation temperature at the fin apex, 
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10. zero film thickness at the fin apex, and 
ll. negligible curvature of the condenser wall. 

Purnomo [1] developed a two-dimensional Finite Element 
Method solution using a linear triangular finite element 
model as shown in Figure 4. Purnomo modified Corley's assump- 
tion that the fin apex was at the saturation temperature and 
allowed the value of the temperature at the apex to float. 
He assumed a parabolic temperature distribution along the 
mem surface. 

Purnomo's statement of the problem for the formulation 


of the Finite Element Method as shown in Figure 5 is 


a = 0 (6) 


Peet the boundary conditions: 





Ole 2 = 
a) along boundary 1, -k 7 hy (T - ie 
aT _ 
ie «along boundary 3, =| ae hy (T - et) 
; ae 
¢) along boundaries 2 and 4&, ae 0 


A detailed description of the numerical formulation is pre- 


Senced in his thesis. 
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Purnomo's computer program consisted of a main program 
mma three subroutines; 

a) the routine "COORD" need to define positions of 
system coordinate points, 

b) the routine "FORMAF" used to formulate the 
Finite Element Method equations, and 

c) the routine "BANDEC" as an equation solver for a 
symmetric matrix which has been transformed into 


banded form. 


Seer, COMPUTER PROGRAM DEVELOPMENT 

Purnomo's [1] two-dimensional finite element program 
is the basis for the present analysis code. The first 
task undertaken in the development of this thesis was to 
check Purnomo's code for validity. Theoretical results 
using Purnomo's code for the case of zero fins were compared 
to theoretical results and experimental data for a smooth 
tube obtained by Schafer [6]. Purnomo's results exceeded 
Schafer's results by a factor of approximately two. In 
studying this discrepancy, an error was discovered in 
Purnomo's code. In encoding formula (II-11) [1] for mass 
Flow rate in FORTRAN, the' sin @ ' term was dropped. The 
correct form of the equation is shown below. 

ow (Ry + xsin eos Sin ¢ 


ee (x) = ———_+-—- *(x)etd #(X) “tan a) 
3 Us 


Cc 
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The effect of this error was that the film thickness along 
the length of the condenser remained very small. Purnomo's 
original code consequently indicated a much higher value of 
heat transfer rate than the correct analysis did. Ina 
correct analysis code, the condensate film thickness grows 
continuously until the condensate reaches the evaporator. A 
plot of the condensate film thickness calculations from 
Purnomo's original code for the case of a smooth tube, his 
corrected code, and Schafer's program for a smooth tube are 
mew in Figure 6. Purnomo's corrected code for the case of 
zero fins agreed to within 8 percent of the results for 

heat transfer rate using Schafer's [6] analysis for a smooth 
mabe . 

Since no experimental data existed for further comparison 
of the finned model, the next task undertaken was to adapt 
the analysis code to permit automated design and sensitivity 
analysis using COPES/CONMIN. Many modifications were made, 
some of which are mentioned here. The program was rewritten 
in subroutine form and segmented into input, execution, and 
output sections to make it compatible with COPES. Since 
COPES was written to use single precision mathematics and 
the analysis code uses double precision to allow for possible 
ill conditioning, the subroutine ANALIZ also makes the trans- 
formation from single to double precision. The initial value 


Seeeene tilm thickness is calculated within the code based on 
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a formulation by Sparrow and Gregg [13]. Previously a 
constant value based on the same analysis was used. 

The modified code can be used alone for analysis of a 
given geometry or it can be used in conjunction with COPES/ 
CONMIN for a single analysis, sensitivity analysis, or 
automated design (optimization). 

A listing of the revised computer program is included 


as Appendix B. 


Mee oENSLITIVITY ANALYSIS 

The sensitivity analysis feature of COPES/CONMIN was used 
fjemeobtain data for the heat transfer rate Q as a function 
of various design variables. All design variables were held 
constant except the one of interest which was varied as 
Beeeirtied in block Q of the COPES input data. 

For example, to obtain the heat transfer rate Q for 26 
different values of rotational speed, the following data 
input is required: 

$ Block P 


i 
13 


meBlock Q 
m5, 26 
ee ie, eZ. , 300., 400., S00., 600., 700., 
800., 900., 1000., 2000., 3000., 4000., 5000., 6000., 
etree OmO., L0000.;, L1I000., 12000., 13000., 
mooo0., 18000., 3000. 


ee 





For each plot only the variable on the ordinate is 
changed. All other design variables remain constant. The 
basic design used for sensitivity plots has the design 


variable values shown below: 


eae tL 


Pool ewVno tele soenoltliyitTy ANALYSIS 


condenser length CLI =e oe40 kamenes 
cone half angle CANGL = 1.uU DEGREES 
condenser radius RBASEL = 0.775 inches 
wall thickness Memeriy = 0.03125 inehes 
mim height B ="0-.025 inches 
speed of rotation RPM = 3000 RPM 
Saturation temperature TSAT =e Lr 
ambient temperature TINF Moi SE 
outside heat transfer 2 
@oerricient HINF SeougO btuvHR get) oF 
fin half angle FANGL = 10 
Beati1o Of trough width 
to fin base width ZOA tel 22.6 
number of fins ZFIN = 40.0 


Note: In the code, ZFIN appears on the left hand side of 
Smmecuadtion and is therefore by strict definition not a 
"design variable". ZFIN is calculated from the values of 


@eeoeb, EZERO, and EPSO. 


Pee leolGN OPTIMIZATION 
Counting parameters such as external heat transfer 


coefficient HINF, there are thirteen possible design 
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variables. Nine of these are geometric or functional para- 
meters such as wall thickness, fin height, and speed of 
rotation. The design variables, possible constraint func- 
tions, and the objective function appear in the Global 

common block GLOB1 in the code and are listed below by fortran 
name for clarity. 


Reo ecN VARTABLES 
er 
CANGL 
ferro 
Ree 
Toes 
BFIN 
TZ 
TSS 
NE 
Pet © 
FANGL 
ZOA 
REM 


CONSTRAINT FUNCTIONS 
ZFIN 
BOA 
Darn 


OBJECTIVE FUNCTION 
QTOT 


There are a wide variety of design problems that can be 


pursued with the code. For example one might wish to 
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determine the smallest length condenser section and the best 
internal geometry for a specified heat transfer rate-- 
perhaps to cool an electric motor. 

To maximize heat transfer rate through the condenser 
wall the designer adds a number of fins to increase the 
inside surface area. As more and more fins are added 
however, the cross-sectional area for conduction through 
each fin is decreased. Also, the film thickness of the 
condensate in the though increases and in fact could com- 
pletely cover the fins and substantially reduce heat trans- 
fer through the fin. So there should exist some optimum 
combination of number of fins, fin height, fin half angle, 
Sm@a@mratlo Of trough width to fin width that will permit 
maximum heat transfer rate. 

The design study undertaken was to determine the fin 
Metent, fin half angle, and fin spacing which would yield 
the maximum heat transfer possible. The design variables 
then were BFIN, FANGL, and ZOA. Other potential design 
variables were held constant. The objective function to be 
maximized was QTOT, heat transfer rate out of the condenser. 

For comparison, the theoretical upper limit on heat 
transfer was calculated based on an external surface tempera- 
ture equal to the working fluid saturation temperature. This 
assumes that there is no thermal resistance across the 
condensate and the condenser wall. When this upper limit 


was used, the maximum heat transfer rate was predicted to be 


SS 





memeo22 BrU/HR using the following formula: 


where 


ay 


L 


Max 


=h 2mr L (Twall - To) (8) 


Sines Gemcemvectmuve heat transrer coefficient 


2 oP) 


Cso@o, SL U/HR- Fr 
average outside radius of condenser wall (0.0056 ft) 


condenser length (0.75 feet) 


Twall=temperature of the outside wall (100°F) 


Too 


ambient temperature (60°F) 


Certain constraints were placed on the design based on 


engineering judgement. For example, the number of fins was 


not allowed to exceed 400 and the minimum fin half angle 


allowed was 10 degrees. These values were based on structural 


and manufacturing considerations. 
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Ve ecco brs 


Pee LN TRODUCTION 

myencomputer code was used in conjunction with the COPES/ 
CONMIN program for sensitivity analysis and design optimiza- 
tion for maximum heat transfer rate. Numerical results are 


discussed below. 


Soe NSLITIVITY ANALYSIS 

As shown in Figures 7 and 8, the heat transfer rate in- 
creases for an independent increase of the design variable 
of interest, but levels off at the theoretical maximum heat 
transfer rate for the given overall geometry. The heat trans- 
fer rate then is limited by the external resistance which has 
Meecome the controlling factor. 

In Figure 9, the heat transfer rate increases in a simi- 
lar way for an increase in external heat transfer coefficient. 
Again, the rate of increase appears to lessen for external 


Ee eee 


heat transfer coefficients above 10,000 BTU/HR-FT 
other limiting resistances. In Figure 10 the heat transfer 
rate 1s observed to increase linearly with condenser length. 
This is expected since Q is a function of the area and for small 


values of » the area varies directly with length. In Figure ll, 


Q is seen to rise in a non-linear manner for an increase in 
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cone half angle, ¢. This is because the internal heat trans- 
jer coefficient, hs is a non-linear function of cone half 
angle. 

Purnomo [1] concluded that the heat transfer rate continu- 
ously increased as the fin half angle decreased and that this 
was largely a result of the fact that the number of fins in-~ 
creased at the same time. He also stated that the increase 
in heat transfer was only slight when the fin half angle was 
less than 11 degrees. This iS somewhat misleading in that in 
his analysis the number of fins was being changed with every 
change in fin half angle. This author has drawn a different 
conclusion. Figure 12 shows a plot of heat transfer rate vs. 
varying fin half angle all for a condenser with 40 fins. The 
heat transfer rate as a function of fin half angle rises 
sharply from 1-11 degrees and continues to rise, but less 
steeply, as the fin half angle increases beyond 11 degrees. 
ieemmaximum fin half angle possible with 40 fins is 67.5 


degrees. 


Sec ONSTRAINED OPTIMIZATION 

In the design problem undertaken to determine the optimum 
internal geometry for maximum heat transfer, numerous runs 
were made for condensers made of copper, stainless steel, and 
a ceramic material. These materials have thermal conduc- 
[very Values of 231,9, and 1.0 BIU/HR-FT-°F respectively. 


It was expected that for each material a different optimum 
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design would emerge. The results, therefore, were unexpected. 
From different starting points (original designs), some out- 
side the feasible region, and for external heat transfer 


noe the same optimum 


coefficients from 1000-50,000 BTU/HR°FT 
design for maximum heat transfer was reached, which is out- 
lined in Table III below. Each material, of course, had a 


different heat transfer rate even though the geometry for 


Maximum heat transfer was the same. 


1S cee et I 


Coo wee D OPTIMIZATION RESULTS FOR COPPER 
Se eileecs Sthal, AND CERAMIC MATERIAL 


Optimum Design for All Materials 


fin height 0.023 inches 
mame half angle 10.0 degrees 
mumoes Of Lins 400.0 

ZotiG Of trough width Os 


tO fin base width 


Nexternal HEAT TRANSFER RATE 

MATERIAL (BTU/HR+ FT? °F) (BTU/HR) 
Copper 1000 12 3807 
50,000 296,560 

Stainless Steel 1000 i 6 08 
50,000 38,650 

Ceramic Material 1000 38uu 
50,000 5237 


og 





meee UNCONSTRAINED OPTIMIZATION 

The results of the constrained optimization runs were so 
unexpected that it was decided to remove the constraint on 
mye number of fins and the trough width to fin width ratio 
and repeat the optimization runs. If another identical design 
was again reached for all materials, the code would have to be 
considered in error. The results are presented in Table IV. 

The optimum design for maximum heat transfer was different 
for each material and each heat transfer coefficient. For 
each material, the optimum fin height b increases as the ex- 
ternal heat transfer coefficient increases. This provides 
more fin surface area to increase heat transfer rate as the 
Outside heat transfer resistance no longer controls. Despite 
a fifty-fold change in heat transfer coefficient, each 
material maintained essentially a constant 6*/b ratio. The 
higher conductivity materials required less exposed fin sur- 
face than the lower conductivity materials. Copper, for 
instance, could be 57 percent covered by the trough condensate 
while the ceramic material was only 14-19 percent covered for 
maximum heat transfer. 

In each material the number of fins decreased for an 
increase in external heat transfer coefficient. This provides 
more space in the troughs to carry the increased condensate 


which results. Otherwise, fin performance would be degraded. 
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For an increase in external heat transfer coefficient the 
highest percentage change in heat transfer rate occurs in the 
material with the highest thermal conductivity. In the 
ceramic material, for instance, a fifty-fold increase in ex- 
ternal heat transfer coefficient results in only a 36 percent 
increase in heat transfer rate. In the copper condenser, the 
heat transfer rate 1S increased by a factor of 21. This shows 
the strong dominance of the wall resistance to heat transfer 
in the ceramic material. 

For all three materials, the trough 1s essentially eli- 
minated for the optimum design. Why the trough 1s not 
entirely eliminated in the stainless steel condenser section 
me not clear. 

Because the manufacture of a very large number of fins in 
a small diameter condenser is not practical, it is desirable 
to use the constrained optimization results as a design guide- 
line. 

A comparison of the unconstrained and constrained design 
for a copper condenser section, for instance, shows that 
using the more realistic number of 400 fins and its proper 
fin height b instead of 2040 fins results in a heat transfer 
rate difference of only 0.6 percent, This is shown clearly 


fie lable V. 
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TABLE V 


Comparison of Unconstrained and Constrained Optimum Design for 
a Copper Condenser Section with External,Heat Transfer Coeffi- 


Gtemeronresooood BLU/HR-FT -°F. 
Number Difference 
b Cinches) of Fins Ome lO, ak > in Q 
constrained O02 3 400 Ze. 560 ---- 
unconstrained OiG19 2040 Zoom so 0 0.6% 


The author's conclusion is that within the realm of what 
can now be manufactured, the same basic design is best for 
all materials regardless of the external heat transfer 
coefficient. That is, the designer should machine as many 
Mmemeeas the condenser material and the manufacturing process 
will allow. The code should be used to determine the fin 
height to avoid degrading the fin efficiency by too high a 
level of condensate in the trough. As seen in Table V, it 
may be possible to lower the constraint on the number of fins 
more than once and compare the resulting heat transfer rate 
to achieve an effective but less expensive design to manu- 


facture. 


fee CAUTIONARY NOTE 

The reader should be cautioned that when an analysis code 
based on assumptions made for certain conditions is linked to 
an optimizer, the code may change the geometry or other 


conditions such that the original assumptions are no longer 
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valid. Judicious use of side constraints and engineering 
common sense can prevent disaster. Numerical solutions must 
Hee be blindly accepted. 

Merscnould be noted that all calculations in this thesis 
were made with water as the working fluid. Sensitivity data 
were run with an original rotational speed of 3000 RPM while 
design optimization was done for a rotational speed of 3600 
RPM. Different’combinations of operating parameters, working 
fluids, and condenser section materials may not behave in a 
predictable manner. While certain trends are shown for the 
particular condenser section studied here, it would be 
preferable to modify the code to analyze the particular 
problem at hand than to extrapolate these findings to tao 
broad an application. 

While a design improvement is almost certain using an 
optimization routine like COPES/CONMIN, there is no guarantee 
that a global optimum has been found. Therefore, for a given 
design problem, the designer should start the optimization 
process from several different original designs and compare 


the results. 
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V. CONCLUSIONS 


1. The computer code is valid and can be used for single 
analysis, sensitivity analysis, and automated design of an 
internally finned rotating heat pipe. It can also be used to 
generate numerical data to find a correlation useful for 
design. 

feet or zero fins, the code converges to the results for 
a smooth tube obtained by Schafer [6]. 

3. The heat transfer rate increases for an independent 
increase in fin half angle, rotational speed, and number of 
fins but levels off at the theoretical maximum heat transfer 
rate for the heat pipe. The heat transfer rate continues to 
increase for an increase in condenser length and cone half 
angle. 

u. Within the realm of known materials, maximum heat 
transfer occurs for the same fin geometry regardless of the 
external heat transfer coefficient. That 1s, for a given 
condenser radius, one should machine as many fins as possible 
with the condenser material used to maximize heat transfer. 
The code should be used to determine the correct fin height 
to avoid degradation of fin performance by too high a level 


of condensate in the trough. 
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5. When the constraint on minimum fin half angle is re- 
moved the optimum design is a large number of very thin fins. 
A practical design then might be a series of long rectangular 


mans. 
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VI. RECOMMENDATIONS 


Based upon the calculated results of this thesis, the 


following recommendations are made. 


i 


Build an internally: finned heat pipe to obtain 
experimental data for comparison to the analytical 
predictions. 

Generate a new model to analyze the case where cone 
half angle is zero. This would greatly decrease the 
cost of manufacture, and permit the use of internally 
finned heat exchanger tubing which is commercially 
available. 

Analyze different shaped fins, including rectangular. 
Obtain numerical data using the code to attempt to 
find a correlation for the heat transfer rate as a 
function of important variables: 

Q = Q (temperature, condenser geometry, fluid 
properties, condenser material properties, 
heat transfer coefficients, etc.) 

Modify the code to include as constraint functions 
structural failure modes such as rupture due to the 


internal spinning mass, buckling, and whirling. 
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Figure 1. Schematic Drawing of a Rotating Heat Pipe 
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meeure 2. Operating Limits of a Typical, Water-Filled 
Rotating Heat Pipe 
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Figure 3. Internally Finned Condenser Geometry, 
Showing Fins, Troughs, and Lines of 
symmetry 
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Figure 4. Condenser Geometry Considered with 25 
Linear Triangular Finite Elements 
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Mesure 5. Differential Equation and Boundary Conditions 


Considered in the Analysis of Purnomo [1] 
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Figure 7. Heat Transfer Rate (Q) vs. Number of Fins 
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Fagure 8. Heat Transfer Rate (0) vs. RPM 
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Figure 9. Heat Transfer Rate (Q) vs. External Heat 
Transfer Coefficient 
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Figure 10. Heat Transfer Rate (Q) vs. Condenser 
Length 
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Figure ll. Heat Transfer Rate (Q) vs. Cone Half Angle 
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APPENDIX A 


USER'S MANUAL 


This ees x describes the data cards for the use of 
mae, computer program. 

The data is divided into the COPES/CONMIN section and the 
heat pipe analysis program section. 

The COPES data 1s segmented into "blocks" for convenience. 
All formats are alphanumeric for TITLE and END cards, F10 for 
real data, and I10 for integer data. Comment cards may be 
inserted anywhere in the data deck prior to the END card and 
are identified by a dollar sign ($) in column one. The COPES 
data deck must terminate with an end card containing the word 
ane in columns 1-3. 

Information is included in this appendix pertaining to 
data needed for single analysis, sensitivity analysis, or 


Metimazation using COPES/CONMIN. 


UNFORMATTED DATA INPUT 

While the user's sheet defines COPES data in formatted 
fields of ten, the data may actually be read in a simplified 
fashion by separating data by commas or one or more blanks. 
If more than one number is contained on an unformatted data 
card, a comma must appear somewhere on the card. If exponen- 
tial numbers such as 2.5+10 are read on an unformatted card, 
there must be no embedded blanks. Unformatted cards may be 


intermingled with formatted cards. Real numbers on an 
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unformatted card must have a decimal point. 


EXAMPLES 


Unformatted data; 

oye. 3,).0+20,-5.1 

fees, 1 .0+20,,-5.1 

meds 1.0+20,, -5.1 

5 Pees, Le UtZ20 0 =o .1 

Equivalent formatted data; 

Sole 10 20 30 40 20 60 70 


. 7 ye eee Q -5.1 
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DATA BLOCK A 


MmeerRiPT LON: COPES TITLE CARD 


FORMAT : 20A4 
Ee Zz 3 + S 6 7 8 
Title 
PLELD CONTENTS 
1-8 Any 80 Character title may be given on this card. 


GZ 





DATA BLOCK B 
DESCRIPTION: 


FORMAT : 


BLELD 


PAEREO 


PROGRAM CONTROL PARAMETERS 
Z 3 uy 5 6 7 8 
NDV NSV 
CONEENTS 
NCAP Cateculation control 
O- Read input and stop. Data of blocks A,B 
and V is required. Remaining data is 
Oprlona ls 
l- One cycle through program. The same as 
executing ANALIZ stand-alone. Data of 
blocks A,B and V is required, remaining 
data is optional. 
7 Oper caetoma —Dacaros, Blocks of A-~I and V 
is required. Remaining data is optional. 
3- Sensitivity analysis. Data of blocks, 
A,B,P,Q and V 1s required. Remaining data 


SiG ieee atl. 


NDV: Number of independent design variables in 


NSW: 


Opeimaza elo. 
Number of variables on which sensitivity 
analysis will be performed. 
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DATA BLOCK C 
DESCRIPTION: 


mORMAT : 


BLELD 





Omit if NDV=0 in block B 


INTEGER OPTIMIZATION CONTROL PARAMETERS 


PVE OSI CNDIR ee NSsCSh ITRM LINOBJ NACMXI 


CONTENTS 


tl wit GCOnerOoL used in the optimization 
program CONMIN. 

f-wiiememimMtedurinca Optimization. 

eevee eainad final Optimization 
information. 

2- Print above plus objective function value 
and design variable values at each iteration. 

3- Print above plus constraint values, direction 
vector and move parameter at each iteration. 

4- Print above plus gradient information. 

5- Print above plus each proposed design vector, 
objective function and constraint values 
during the one dimensional search. 

ITMAX: Maximum number of optimization inter- 
actions allowed. Default=20. 

ICNDIR: Conjugate direction restart parameter. 
GT. O- Scale design variables to order of 
magnitude one every NSCAL iterations. LT. 0- 
Scale design variables according to user in- 
put scaling values. If not zero, NDV + 1 
1s recommended. 

ITRM: Number of consecutive iterations which 
must satisfy relative or absolute convergence 
criterion before optimization process is 
terminated. Default=3. 

LINOBJ: Linear objective function identifier. 
If the optimization objective is known to be 
a linear function of the design variables, 
set LINOBJ=1. Default=Non-linear. 
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NACMX1: One plus the maximum number of active 
constraints anticipated. Default-NDV=2. 
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DATA BLOCK D Omit 1f NDV=0 in block B 
BeecktPTiON: FLOATING POINT OPTIMIZATION PROGRAM PARAMETERS 


BORMAT: Ee 


mC H FUCHM ear Crain Cie LUMEN ~ THETA 


NOTE: Two cards are read here. 


mot LD CONTENTS 





He FDCH: Relative change in design variables in 
calculating finite difference gradients. 
DetauiG= 0.01 - 


2 Poet Mimi atosOlutemstep in finite dif- 
ference gradient calculations. Default= 
Ceo 
3 CT: Constraint thickness parameter. Default= 
=-0.05. . 
4 CTMIN: Minimum absolute value of CT consi- 


dered in the optimization process. 
Default=0.004. 


2 CTL: Constraint thickness parameter for 
Pimcameeonormaamers., » Default=-0.01. 
6 CM ss iianimum aosolute value of CTL con- 


sidered in the optimization process. 
beraunit=0 001. 


7 THETA: Mean value of push-off factor in the 
method of feasible directions. Default= 
oO. 
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DATA BLOCK D mee lf NPV=0 in biock B 


BORMAT : 4FLO 


DELFUN DABFUN ALPHAX ABOBJ1 


ELELD CONTENTS 


1 DELFUN: Minimum relative change in objective 
function to indicate convergence of the 
optimization process. DEFAULT=0.001l. 

2 DABFUN: Minimum absolute change in objective 
function to indicate convergence of the 
OPulmuzidulomupeceecess.. DEFAULT=0.00ae times 
the initial objective value. 

3 ALPHAX: Maximum fractional change in any design 
variable for first estimate of the step in 
the one-dimensional search. DEFAULT=0.1. 

u ABOBJ1: Expected fractional change in the 
objective function for first estimate of the 
step in the one-dimensional search. 
DEFAULT=0.1. 


REMARKS 


1) The DEFAULT values for these parameters usually work well. 
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DATA 


BmOcK -E Onmeeeir NDV=0 in block B 


MmeockEPTION: TOTAL NUMBER OF DESIGN VARIABLES, DESIGN 


Sreeel lye CDENTITICATION AND SIGN 


monMAT ; Zeon E10 


fee LD 


NDVTOT IOBJ SGNOPT 


CONTENTS 


NDVTOT: Total number of variables linked to the 
design variables. This option allows two or 
more parameters to be assigned to a single 
design variable. The value of each para- 
meter is the value of the design variable 
times a multiplier, which may be different 
for each parameter. DEFAULT=NDV. 

IOBJ: Global variable location associated with 
the objective function in optimization 

SGNOPT: Sign used to identify whether function 
PSEC Demat Zed,or Manimized. +1.0 indi- 
cates maximization. -1.0 indicates minimi- 
zation. If SGNOPT is not unity in magnitude, 
it acts as a multiplier as well, to scale the 
magnitude of the objective. 
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Bara BLOCK F 


Beach lI PIrION: DESIGN VARIABLE BOUNDS, 


SCALING FACTORS 


BORMAT : 4FLO 


VLB VUB x SCAL 


ieee: MOV=0 an block B 


INITIAL VALUES AND 


See READ ONE CARD FOR EACH OF THE NDV INDEPENDENT DESIGN 


VARIABLES. 

pees LD CONTENTS 

ik VLB: Lower bound on the deSign variable. If 
low tee or flo weno. 1Owem: bound. 

2 VUB: Upper bound on the design variable. If 
Jeno mire. NO .wpper bound. 

3 X: Initial value of the design variable. If 
X is non-zero, this will supercede the value 
initialized by the user-supplied subroutine 
ANALIZ. 

4 SCAL: Design variable scale factor. Not used 


Pe oe Momo tile hlOCk —C. 
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mara BLOCK G Omit if NDV=0 in block B 
DESCRIPTION: DESIGN VARIABLE IDENTIFICATION 


Meme: §=2L10,F10 


NDSGN IDSGN AMULT 


fetes READ ONE CARD FOR EACH OF THE NDVTOT DESIGN VARIABLES 
Rielnh offbeOnDER AS IN BLOCK F. 





ELELD CONTENTS 

a NDSGN: Design variable number associated with 
this variable. 

2 IDSGN: Global variable number associated with 
this variable. 

3 Livi cOtctanrtemiltiplier on this variable. 


ion 7 Moment etnoe variable will be the value 
of the design variable, NDSGN, times AMULT. 
DEEAULI=Hal.0. 
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DATA BLOCK H Titel E=NbY 0 in block 8B 
meee RI PI LION: NUMBER OF CONSTRAINED Peeve te eens 


SeevwAT: I10 


al: 2 3 4 5 6 7 8 
NCONS 
be LD CONTENTS 
1 NeEOiec » Number Of conStraint sets in the 


optimization problem. 


REMARKS : 
1) If two or more adjacent parameters in the global common 


block have the same limits imposed, these are part of the 
same constraint set. 
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DATA BLOCK: [ eeeeeiov=O0 an block B;, or NCONS=0 in 
Dock 


MeockRLPTION: CONSTRAINT IDENTIFICATION AND CONSTRAINT BOUNDS 
monMAl; 31f10 


aCON JCON LCON 


MOTE: READ TWO CARDS FOR EACH OF THE NCONS CONSTRAINT SETS. 


Bee LD CONTENTS 

i ICON: Farst global number corresponding to the 
Gonstraint Set. 

2 JCON: Last global number corresponding to the 
Sonstraimr set. DEFAULT =ICON. 

3 mOONemEtrnear CoOnstrartnt 2dentiftier for this 
constraint set. LCON=1 indicates linear 
constraints. 
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meta BLOCK: I COON TEIGED) 


FORMAT: 4F10 


BL SCALIL BU SCAL2 


foe LD CONTENTS 


i. BL: Lower bound on the constrained variables. 
Pao ee — leet lS eeno lower bound. 

2 SCAL1: Normalization factor on lower bound. 
DEFAULT=MAX of ABS(BL), 0.1. 

3 BU: Upper bound on the constrained variables. 
Beep t sot Wer lo sno upper bound. 

4 SCAL2: Normalization factor on upper bound. 


DEFAULT=MAX of ABS(BU), O.1. 


REMARKS : 


1) The normalization factor should usually be defaulted. 

meeerne Constraint functions sent to CONMIN are of the forn: 
Meme vAbUE)/SCALL -LE. 0.0 and (VALUE - BU)/SCAL2 .LE. 
Oe 0 . 

3) Each constrained parameter is converted to two con- 
constraints in CONMIN unless ABS(BL) or ABS(BU) exceeds 
[ere ris, in whieh cease no constraint is created for that 
meund. 
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merA BLOCK: P Gime eretov—-0 in block B 
MescRIPTION: SENSITIVITY OBJECTIVES 


MenMAT: 2110 


m20BJ IPSENS 


Mere LWO OR MORE CARDS ARE READ HERE. 


FIELD CONTENTS 





i NSOBJ: Number of separate objective functions 
to be calculated as function of the sensi- 
tivity variables. 

2 ieomiio. tenes control. [fePPSENS=GT.0, de= 
tailed print will be called at each step in 
the sensitivity analysis. DEFAULT=No 
ete de 
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BATA BLOCK: P (CONTINUED) 
DESCRIPTION: 


MoRMAT: 8110 


NSN1 NSN2 NSN3 NSN4Y 


ELELD CONTENTS 


1-8 NSNI: Global variable number associated with 
the sensitivity objective functions. 


REMARKS: 


1) More than eight sensitivity objectives are allowed. 
data cards as required to contain data. 


ess 





DATA BLOCK: Q ite nh Noy =0 in block B 
MeocREPLION: SENSITIVITY VARIABLES 


Beet: 2110 


foeN> NSENS 


Meme «60CREAD ONE SET OF DATA FOR EACH OF THE NSV SENSITIVITY 
Perea epiboew .JWO OR MORE CARDS ARE READ FOR EACH SET 


OE DATA. 
mere LD CONTENTS 
i ISENS: Global variable number associated with 
the sensitivity variable. 
P NSENS: Number of values of this sensitivity 


variable to be read on the next card. 
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BatA BLOCK: Q (CONTINUED ) 
MmeocREPTION: 8sF10 


PORMAT : 


SNSL SNS 2 SNS 3 SNS4Y 





BLE LD CONTENTS 

1-8 SENSI: Values of the sensitivity variable. 
PT episaio, sli correspendas to the nominal 
value. 

REMARKS : 


1) More than eight values of the sensitivity variable are 
allowed. Add data cards as required to contain the data. 
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DATA BLOCK: V 
DESCRIPTION: COPES DATA 'END' CARD 


monMaAr: 3Al 


co 2 3 Ly 5 6 7 8 
END 
PLE LD CONTENTS 
i. muceword “END san columns 1-3. 


REMARKS 
ieeetnis card MUST appear at the end of the COPES data. 


Mmmeenis ends the COPES input data. 
3) Data for the user-supplied routine, ANALIZ, follows this. 


ES 


Seem PIPE ANALYSIS 
Data for the heat pipe analysis follows the 'END' card 

in the COPES data deck. If the general design capability of 
COPES/CONMIN is not needed, the heat pipe analysis can 
be run by setting NCALC = 1 in field number 1 of data block 
B; Or in a stand-alone mode by using the following main 
program. 

eerie PROGRAM FOR HEAT PIPE ANALYSIS 

Cee eee crncULe, AND PRINT RESULTS 

Dew hipere Ahoy = 1.3 
mer CALL ANALIZ CICALC) 


S20r 
END 
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DATA BLOCK: AA 
DESCRIPTION: ELEMENT CONNECTIVITY 


BORAT: 315 


NEL NSNP NBAN 


PLE LD CONTENTS | 

il NEL: Number of elements 

2 NSNP: Number of system nodal points 
3 NBAN: System band width 
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DATA BLOCK: BB 


mescnktPTION: ELEMENT CONNECTIVITY 


FORMAT: 415 eo ee =e NEL 


mei LCON@IEL,1) 


met 1.D CONTENTS 
i IEL: The element number 
2 Tee lela): “System modal point corresponing 


to nodal point I of element TEL 


REMARKS : 

1) Number all elements with convective boundaries first from 
top to bottom, then number the remaining elements. 

2) Number the nodal points of each element moving in a 
counter-clockwise direction. 

3) The elements with convective boundaries have nodal points 


l and 2 located on the convective boundary. 
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DATA BLOCK: CC 
DESCRIPTION: CONDENSER GEOMETRY 


MeeMAT: /G10.5 


Gist CANGL RBASEL ioe «LECT 6c BF IN az 


BLELD CONTENTS 


CLI: Condenser length (inches) 

CANGL: Cone half angle (degrees) 

RBASEI: Inside radius of condenser base 
(inches) 

R2I: Intermediate radius (inches) 

THICKI: Condenser wall thickness (inches) 

berN=  Hedgehe of fin Cinekes) 

- TZ: Nodal point temperature initial guess 

(degrees F) 


Wh / 


“IO on 


REMARKS 


J) Set TZ equal to TSS-Saturation temperature of the 
working fluid 


8 2 





DATA BLOCK: 
mesCRIPTION: 


FORMAT: 


PLELD 


Pevi lee le MENT GEOMETRY 


Nao. NEFB NBOTI NBOTF 


CONlTEN iS 


NDIV: Number of increments along the length 
of the condenser. 
NEST: Number of the element on the right end 
See el tet On ie 
NEFB: The element number with convective 
boundary located at the base of the fin. 
NBOTI: The element number with convective 
beundary Located at the right hand of the 
bottom side. 
NBOTF: The element number with convective 
boundary located at the left hand of the 
bottom side. 
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DATA BLOCK: EE 
DESCRIPTION: DATA FOR RUNNING 


BerMAT: 4F10.2 


REM TSS eat HINF 


BLELD CONTENTS 

1 RPM: Rotation rate of heat pipe (RPM) 

2 TSS: Saturation temperature of the working 
Fluid (degrees F) 

3 TINF: Outside temperature (degrees F) 

4 Hive Oueside convective heat transfer 


Goctrterenr (SlU/HREY °F) 
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DATA BLOCK: FF 
DESCRIPTION: CONVERGENCE CRITERION 


BeBMAT: G10.9 


HE 2 3 4 5 6 7 8 
cet 
Ele LD CONTENTS 
iL CRIT: Convergence criterion on finite element 


SoluEtoOnwt ore imeremental Neat transter rate 
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DATA BLOCK: GG 
DESCRIPTION: INTERNAL FIN GEOMETRY 


FORMAT: 2G10.5 


1 2 3 4 2 6 di 8 
FANGL ZOA 
fp LELD CONTENTS 


FANGL: Fin half angle (degrees) 
ZOA: Ratio of trough width to fin base width 


Nr 
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DATA BLOCK: HH 


MescRLPTION: INTERNAL FIN GEOMETRY 


moRMAT: 15 
1 "2 3 4 5 6 7 8 
ine 
pele 1D CONTENTS 
il IFF: (n-1), where n is the number of rows of 


the upper triangular fin section 


NOTE: See Figure A-1. 
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DATA BLOCK: II 
DESCRIPTION: INTERNAL FIN GEOMETRY 


FORMAT: i615 b=, 1FF 


MEENCI) KEECT) 


FIELD CONTENTS 


1 KFIN: The number of system nodal points located 
on the symmetric boundary of triangular fin 
section, but does not include the system 
nodal points located at the base of the fin 
and the apex. 

2 KFF: The number of system nodal points located 
along the fin convective boundary, but does 
not include the system nodal points located 
at the base of the fin and the apex. 


MOTE: See Figure A-1l. 
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DATA BLOCK: JJ 
DESCRIPTION: INTERNAL FIN GEOMETRY 


PORMAT: 2610.5,415 


DOBF DOTH ELE JLC JINT For 





FIELD COTES 

i POlEes  umper of column within fin. 

Z DOTH: Number of column within though. 

3 ieee (i hemmnumoer Of the system nodal point 


located at the junction of the symmetry 
boundary and the line of intersection 
between the fin and the condenser wall. 


4 JLC: The number of the system nodal point 
located at the center of system coordinates. 
5 JINT: The numerical difference between the two 


adjacent system nodal points vertically 
at the condenser section. 
6 aoe new auinber Of Pows WLtenin the wall section. 


MOTE: See Figure A-l, 
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IktF =2 DOBF. = 3.0 ecu =.17 
KFIN=2; 4 DOTH = 1.0 JINT = 5 
hey = 3 7° 6 JIC =7 KT = 2 


Figure A-1. Specification for Input Data to Determine 
the Coordinates of the System Nodal Points 
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